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ABSTRACT 

Molecular line observations of NH 3 (J,K)=(1,1), (2,2) and CH 3 OH at 24.93GHz taken with the Aus- 
tralian Telescope Compact Array (ATCA) toward the massive twin cores NGC6334I & I(N) reveal 
significant variations in the line emission between the two massive cores. The UCHn region/hot core 
NGC6334I exhibits strong thermal NH 3 and CH 3 OH emission adjacent to the UCHn region and coin- 
cident with two mm continuum peaks observed by Hunter et al. (in prep.). In contrast, we find neither 
compact NH 3 nor thermal CH 3 OH line emission toward NGC6334I(N). There, the NH 3 emission is dis- 
tributed over a broad region (> 1') without a clear peak, and we find Class I CH 3 OH maser emission 
with peak brightness temperatures up to 7000 K. The maser emission peaks appear to be spatially as- 
sociated with the interfaces between the molecular outflows and the ambient dense gas. Peak NH 3 (1,1) 
line brightness temperatures > 70 K in both regions indicate gas temperatures of the same order. NH 3 
emission is also detected toward the outflow in NGC6334I resulting in an estimated rotational temper- 
ature of Trot ~ 19 K. Furthermore, we observe CH 3 OH and NH 3 absorption toward the UCHn region, 
the velocity structure is consistent with expanding molecular gas around the UCHn region. Thermal 
and kinematic effects possibly imposed from the UCHn region on the molecular core are also discussed. 
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1. INTRODUCTION 

The prominent molecular cloud/Hn region complex 
NGC6334 at a distance of 1.7 kpc is composed of a series 
of luminous star forming regions at various stages of evolu- 
tion (Neckel 1978; Straw & Hyland 1989). During the ini- 
tial radio and infrared studies completed over twenty years 
ago, a confusing family of nomenclature was introduced 
for this region. Letters A-F refer to VLA 6cm continuum 
sources (Rodriguez et al. 1982), while Roman numerals I 
through VI refer to far-infrared sources (McBreen et al. 
1979; Loughran et al. 1986). While each center of activity 
contains CO hot spots and various masers, the youngest re- 
gion appears to be radio source "F" (which corresponds to 
infrared source "I" and IRAS source 17175-3544) located 
at the northeastern end of the complex. Even the earliest 
far-infrared maps (Gezari 1982) resolved the emission of 
source I into two cores: the sources I and I(N), which are 
separated by ~ 90". With the cold dust emission, only a 
few red infrared sources and weak cm continuum emission 
(e.g., Gezari 1982; Tapia et al. 1996; Sandell 2000; Car- 
ral et al. 2002), Source I(N) is considered to be the least 
evolved object in the region. 



1.1. NGC6334I 

A well-known, cometary compact Hn region lies at the 
heart of source I (e.g., Carral et al. 1997) and contains 
associated OH masers (Brooks & Whiteoak 2001; Gaume 
& Mutel 1987), Class II CH 3 OH masers (Ellingsen et al. 



1996) and H 2 masers (Moran & Rodriguez 1980; Forster 
& Caswell 1989). Recent mid-infrared imaging with OS- 
CIR at the Keck II telescope has identified a candidate star 
(IRS IE) for the excitation source of the compact Hn re- 
gion (De Buizer et al. 2002). Single dish bolometer maps 
from JCMT reveal a compact (< 20") strong (sub)mm 
continuum source with an estimated dust temperature of 
100 K (Sandell 2000). A SEST line-survey (Thorwirth 
et al. 2003) shows rich molecular line emission in the mm 
wavelength range similar to prototypical high-mass star- 
forming regions such as Orion-KL. Comparable results are 
obtained in a mini line survey near 345 GHz (McCutcheon 
et al. 2000). A CO(2-l) outflow with a velocity range 
of - 150 km s^ 1 emanates from NGC6334I (Bachiller & 
Cernicharo 1990) and a similar outflow is also seen in the 
H76a and H92a recombination lines (de Pree et al. 1995). 
Several knots of shock-excited H 2 emission have been de- 
tected near the ends of this outflow (Davis & Eisloeffel 
1995; Persi et al. 1996; Megeath & Tieftrunk 1999) along 
with NH 3 (3,3) masers (Kraemer & Jackson 1995). The 
driving source(s) for the outflow were claimed to be either 
the powering source of the UCHn region, or an adjacent 
object. One candidate is the cold 20/zm source LKDJ2 
located 6.5" to the north-west of the UCHn region (as re- 
solved by MIRAC2 at IRTF), which is not seen at shorter 
wavelengths (Kraemer et al. 1999). Alternatively, the pow- 
ering source of the UCHn region and LKDJ2 may both 
drive independent outflows. 
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1.2. NGC6334I(N) 

Similar to source I, the JCMT bolometer maps show 
I(N) to be a strong source of (sub)mm continuum emission 
(Sandcll 2000). The estimated dust temperature is lower 
around 30 K. While source I is stronger than I(N) in CS (5- 
4) (Kraemer & Jackson 1999) and in CS (7-6), the HC 3 N 
(15-14) line is much stronger in I(N) than in I (Megeath 
& Ticftrunk 1999; Sollins & Megeath 2004). Further dif- 
ferences can be seen in the 230 GHz line survey where 
the CH 3 OH bands (at 240-242 GHz and 250-252 GHz) are 
fainter in I(N) than I (Thorwirth et al. 2003), probably 
indicating a younger age for this region. Similar conclu- 
sions were drawn by McCutcheon et al. (2000) based on 
their smaller line survey and spectral line mapping of the 
two regions. Indeed, in striking contrast to I, the I(N) re- 
gion shows no mid-infrared emission and only a few red 
infrared sources (Tapia et al. 1996). Traditionally, I(N) 
has also been known as lacking cm continuum emission; 
however, two faint (~ 0.3 mJy), compact cm sources have 
recently been detected with the VLA (Carral et al. 2002). 
One object sits only 7" from the (sub)millimeter contin- 
uum peak, and is consistent with the ionizing radiation 
from a B2 ZAMS star. This cm peak is associated with 
Class II CH3OH maser emission (Ellingscn et al. 1996). 
Furthermore, Class I CH3OH maser emission has been 
mapped at 44 GHz and shows emission peaks distributed 
over approximately 30" (Kogan & Slysh 1998). A bipolar 
outflow has also been mapped in the SiO(2-l) and (5-4) 
transitions with the SEST in south-east north-west direc- 
tion approximately coincident with the distribution of the 
Class I CH3OH maser at 44 GHz (Megeath & Tieftrunk 
1999). These authors also find independent evidence for 
outflow activity in H 2 clumps, which are present but less 
pronounced than in source I. 

1.3. Concerted ATCA and SMA observations 

Because NGC6334I & I(N) are close to each other - thus 
sharing a similar large-scale environment - but in differ- 
ent evolutionary stages, they are the ideal target sources 
to characterize the evolutionary differences in physical and 
chemical properties of massive star-forming regions. Both 
sources have been thoroughly investigated at near/mid- 
infrared wavelength and at lower spatial resolution in the 
radio and mm regime (with some already mentioned ex- 
ceptions). However, high-spatial resolution molecular line 
and continuum studies from short cm to submm wave- 
length have been lacking so far because of its location in 
the southern hemisphere. This situation has changed over 
the last couple of years with the advent of the 12 mm and 
3 mm systems at the Australian Telescope Compact Ar- 
ray (ATCA) and the dedication of the Submillimcter Ar- 
ray (SMA). Therefore, we started a combined effort to in- 
vestigate the molecular line and continuum emission from 
12 mm to the submm regime employing both instruments. 
Here, we present the NH 3 (1,1), (2,2) and the CH 3 OH (near 
25 GHz) results from the ATCA, whereas a separate pa- 
per outlines the first SMA results toward these sources at 
217 GHz (Hunter et al., in prep.). 

The SMA data show three mm continuum sources to- 
ward NGC6334I, two of them adjacent to the UCH11 re- 
gion and one approximately coincident with the cm contin- 
uum peak (Hunter et al. in prep.). Toward NGC6334I(N), 



they detect another mm continuum peak (a few arcsec- 
ond offset from the cm-peak) with a little bit of more ex- 
tended emission. These mm peak positions will be com- 
pared with our NH 3 and CH 3 OH emission throughout this 
paper. Furthermore, Hunter et al. (in prep.) detect very 
collimated SiO(5-4) outflows in NGC6334I & I(N), and 
the outflow orientation in both sources is approximately 
the same in the north-east south-west direction. It is im- 
portant to mention that the new outflow in NGC6334I(N) 
is approximately perpendicular to the previously observed 
larger-scale outflow (Megeath & Tieftrunk 1999). 

2. OBSERVATIONS 

The sources NGC6334I & I(N) were observed on July 
6th and 7th of 2004 with the Australian Telescope Com- 
pact Array (ATCA) in the compact 750D configuration 
with 5 antennas and a baseline range between 31 and 
719 m. Antenna 6 was also included, resulting in a max- 
imum projected baseline of ~ 4.3 km. The phase refer- 
ence centers are given in Table 1. The observed lines are 
the NH3(1,1) and (2,2) inversion lines, with the frequen- 
cies of the main hyperfine components at 23.695 GHz and 
23.723 GHz, and one spectral setup at 24.932 GHz covering 
three CH 3 OH lines at 24.9287 GHz (CH 3 OH(3 2 1 -3i 2 )E), 
24.9335 GHz (CH 3 OH(4 2 , 2 - 4i, 3 )E) and 24.9344 GHz 
(CH 3 OH(2 2 , — 2i i i)E). In addition, we observed the cm 
continuum emission with one spectral unit covering a 
bandwidth of 128 MHz centered at 22.08 GHz. A good uv- 
coverage was obtained through regular switching between 
both sources and the three spectral setups. The spectral 
resolution of the observations was 16 kHz corresponding to 
a velocity resolution of ~ 0.2 km s -1 . The primary beam 
of the ATCA at the given frequency is <~ 130". The data 
were reduced with the MIRIAD package. Applying a ro- 
bust weighting of 1 (closer to natural than uniform weight- 
ing, thus stressing the shorter baselines) the synthesized 
beam is 2.8" x 2.2". We smoothed the data to 0.5 km s -1 
channels, and the rms per 0.5 km s -1 is ~ 6mJy. The 
rms of the 22 GHz continuum and some of the velocity- 
integrated line maps are partly slightly higher, mainly 
because of more pronounced side-lobe effects due to the 
strong emission features. The 3a values are given in each 
Figure separately . 

3. OBSERVATIONAL RESULTS 

The spatial distribution of the molecular gas is signif- 
icantly different between NGC6334I and NGC6334I(N) 
(Figures 1 & 2). Toward NGC6334I, we find strong ther- 
mal, compact NH 3 and CH 3 OH emission in close proximity 
to the UCH11 region. In contrast, the NH 3 emission toward 
NGC6334I(N) exhibits very extended emission without a 
compact component, and we find no thermal CH3OH fea- 
ture but only strong maser emission. 

3.1. NGC6334I 

The NH 3 and CH 3 OH emission from NGC6334I shows 
two emission peaks a few arc-seconds offset from the 
UCH11 region (Fig. 1). The low- level molecular emission 
terminates at the edge of the UCHn region, thus we have 
an example of an observed edge-on orientation between 
an UCHn region next to a hot core. The two molecular 
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line peaks are coincident with two of the three mm con- 
tinuum peaks observed by Hunter et al. (in prep.) with 
the SMA. Interestingly, the molecular line peaks are of 
the same intensity whereas the north-eastern mm-peak has 
about 2.5 times the continuum flux than the south-western 
mm-peak. Another difference between the two peak po- 
sitions is that the north-eastern peak is associated with 
H2O maser emission (Moran & Rodriguez 1980; Forster 
& Caswell 1989) whereas the south-western peak exhibits 
Class II CH3OH maser emission (Ellingscn ct al. 1996). 
More CH3OH, H 2 and OH maser emission is found at 
the interface between the molecular core and the UCH11 
region. Furthermore, we detect extended NH3(1,1) emis- 
sion south-west of the two peak positions. This extended 
emission is significantly weaker in the NH3(2,2) line and 
is barely detected in CH3OH. It is in the same direction 
as the molecular outflow observed in CO by Bachiller & 
Cernicharo (1990) and more recently in SiO by Hunter et 
al. (in prep.). 

The spatial distribution of the NH 3 gas we have ob- 
served varies considerably from the previously reported 
VLA map by Jackson et al. (1988). With a spatial resolu- 
tion of 5.2" x 7.7", they detected more extended emission 
and their peak intensity position is ~ 5" south of ours. 
Tapering our data to the same lower spatial resolution, 
we recover most of their extended structure, but our peak 
positions remains coincident with the two peaks shown in 
Figure 1. Due to the large number of antennas at the 
VLA (27 versus 6 at the ATCA), their uv-coverage should 
be good and allow a reasonable recovery of the large-scale 
structure, but the offset of their peak position remains 
a puzzle. However, our data are less sensitive to large- 
scale emission and therefore trace the small-scale cores far 
better. Additionally, the exact coincidence of our NH 3 
and CH3OH peak positions with the mm peaks reported 
by Hunter ct al. (in prep.) give additional confidence to 
the small-scale gas core distribution presented in Figure 
1. Furthermore, considering the multiple sources we find 
in NGC6334I and the now known outflow in the region 
(Bachiller & Cernicharo 1990), the disk- interpretation of 
the NH 3 data by Jackson et al. (1988) does not hold any- 
more, and it is more likely that their observed large-scale 
emission is somehow associated with the outflow as well. 

3.2. NGC'6334I(N) 

We do not detect any compact thermal NH 3 or CH 3 OH 
emission toward NGC6334I(N) (Fig. 2). However, the 
extended NH 3 emission has about the same peak inten- 
sities as NGC6334I but is distributed over a large area 
(> 1') without a distinct peak position. The mm peak 
detected by Hunter et al. (in prep.) is located approx- 
imately at the center of the large-scale NH3 core. The 
NH 3 (1,1) and (2,2) emission has a similar spatial distri- 
bution as the NH 3 (3,3) emission previously reported by 
Kraemer & Jackson (1999). 

While we do not detect any thermal CH 3 OH emis- 
sion, we find Class I CH 3 OH maser emission from the 
CH 3 OH(3 2 .i - 3i, 2 )E and CH 3 OH(4 2 . 2 - 4i. 3 )E lines 
(nothing from the CH 3 OH(2 2i0 - 2 M )E line). Both 
lines show a maser component at <~ — 3kms~ 1 , but the 
CH 3 OH(4 2 ^ 2 — 4^ 3 )E line exhibits another maser compo- 
nent at <~ — 5kms _1 (Fig. 3). The spatial position of 



the — 3kms _1 components are approximately 7" north- 
east of the mm-peak in the direction of the collimated 
SiO outflow, whereas the — 5kms _1 component is located 
about 10" south of the mm-peak (the exact positions are 
given in Table 2). The observed peak brightness temper- 
ature of the -Skins" 1 line is - 7000 K (Fig. 3). The 
— 5kms _1 component is less clearly associated with one of 
the two outflow axes, but both components are spatially 
close to previously detected Class I CH 3 OH(7o — 61) A + 
44 GHz maser positions (Kogan & Slysh 1998). We note 
that there are significantly fewer, just two, 25 GHz maser 
components detected than 44 GHz components. The 23 
maser spots observed in the latter line are spread over an 
area of 23" x 32", which is somewhat larger, but compa- 
rable to the ~ 15" separation of the 25 GHz components. 
The positional offsets between the 44 GHz and the 25 GHz 
maser positions may be due to the loss of absolute po- 
sitional information during self-calibration of the 44 GHz 
data by Kogan & Slysh (1998). Performing Gaussian fits 
to the CH3OH data in the image plane, we get slightly ex- 
tended features compared to the beam size. This indicates 
unresolved spatial blending of different maser components 
within each channel width of 0.2kms _1 . 

4. ANALYSIS 

4.1. Temperature estimates from the NH 3 observations 

4.1.1. The NGC6334 I peak region 

NH3(1,1) and (2,2) lines are widely used as thermome- 
ters in molecular clouds (e.g., Walmsley & Ungerechts 
1983; Danby et al. 1988). Figure 4 shows NH 3 (1,1) and 
(2,2) spectra toward the two NH 3 peak positions over- 
layed with fits to the hyperfine structure. The hyper- 
fine structure line fits were performed in CLASS (part of 
the GILDAS software package) assuming a Gaussian ve- 
locity distribution, equal excitation temperatures and the 
theoretical intensity ratio and spacing between main and 
satellite hyperfine lines. The fitted parameters are the 
excitation temperature, the opacity of the main hyper- 
fine component, the line- width and the vi sr . Obviously, 
the flat-topped fits do not reproduce the observed spec- 
tra well which is likely due to an extremely high optical 
depth and varying excitation conditions over the width 
of the lines. Missing flux due to missing short spacings 
does not play a significant role since we sample large-scale 
emission in NGC6334I(N) (Fig. 2) and the uv-coverage 
of both observations is nearly the same. These fits do 
not allow us to derive accurate temperatures. Therefore, 
we tried to work only with the satellite hyperfine com- 
ponents of the spectra, and Figure 5 presents the inte- 
grated emission maps of the outer satellite hyperfine com- 
ponents at the most negative velocities (the blend of the 
NH 3 (l,l)F,Fi = 0.5,0-0.5,1 k F,F 1 = 0.5,0-1.5,1 
lines at ~ 23.696 GHz and the NH 3 (2,2)F = 2-1 line at 
~ 23.725 GHz). In the following, we refer to these lines as 
the (outermost) satellite hyperfine (hf) components. An 
intriguing difference between the integrated NH 3 (1,1) and 
(2,2) satellite hyperfine structure line maps is that the 
(1,1) satellite hf-line shows a similar peak intensity dis- 
tribution as the (1,1) main hf-line whereas the (2,2) satel- 
lite hf-line shows an additional intensity peak a few arc- 
seconds south directly at the western edge of the UCHn 
region peak. A (2,2)/(l,l) hf-ratio-map between the two 
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integrated hf-maps stresses this difference as weli (Figure 
5, right panei). The peak NH 3 (2,2)/(1,1) hf-ratios at this 
position arc ~ 4 in contrast to ratios of ~ 0.8 toward the 
two main peak positions. 

The energy levels of the NH 3 (1,1) and (2,2) lines 
are 23 and 64 K above ground, and this line-pair is 
mainly sensitive to temperatures below 50 K. As shown by 
Danby ct al. (1988), rotational temperatures derived from 
NH3(1,1) and (2,2) line observations are not sensitive to 
kinetic temperature changes above about 50 K. Therefore, 
NH3(2,2)/(1,1) hf-ratios > 0.4 are of no use for tempera- 
ture estimates. This is the case for both mm-peak posi- 
tions as well as the additional (2,2) hf-peak close to the 
UCHn region (Fig. 5). However, we observe NH 3 (1,1) peak 
brightness temperatures > 70 K toward the two mm peak 
positions, and these brightness temperatures can be con- 
sidered as lower limits for the kinetic temperatures at these 
positions. As the (2,2) /(1,1) hf-ratio toward the additional 
(2,2) hf-peak is significantly higher (~ 4), the tempera- 
tures there should be higher as well. While it is likely that 
massive protostellar objects at the center of the two mm- 
peaks are causing the high temperatures there, the hot 
spot close to the UCHn region is probably being caused 
by heating of the UCHn region itself. 

The energy levels above ground for the three CH 3 OH 
lines are between 28 and 44 K, thus these lines cannot 
be used for more accurate temperature determinations as 
well. 

4.1.2. The NGC6334 1 outflow region 

The situation is different in the outflow region, and 
we can extract reasonable NH3(1,1) and (2,2) spectra 
toward the NH 3 (1,1) outflow-associated peak at offsets 
(—10", —6"). Fitting the whole hyperfine spectra at this 
position, we can derive a rotational temperature of T rot ~ 
19 K (for details on the temperature calculation from NH3 
observations, see, e.g., Ungerechts et al. 1986). At these 
low temperatures, rotational temperatures correspond well 
to kinetic temperatures (Danby et al. 1988). 

4.1.3. The northern source NGC6334 I(N) 

Since we have no well defined NH 3 peak position in 
NGC6334I(N), we extracted spectra toward the mm peak 
position determined by Hunter et al. (in prep.), which is 
offset by (6.5", 5") from our phase reference center. The 
brightness temperatures are comparable to the brightness 
temperatures in the southern source NGC6334I but the 
opacities are lower as observed in the lower hyperfine line 
strengths (especially in the (2,2) lines, Figure 6). There- 
fore, we can fit these lines reasonably well and again de- 
rive rotational temperatures as done previously for the 
outflow position. The derived NH 3 rotational temper- 
ature is T rot ~ 47 K corresponding to kinetic tempera- 
tures T kin > 100 K (Danby et al. 1988). Thus, again 
we can derive only lower limits for the temperatures to- 
ward NGC6334I(N). This estimate based on the high- 
spatial resolution intcrfcromctric observations is signifi- 
cantly higher than previous estimates of 30 K from single- 
dish NH3 observations with a beam size of 55" (Kuiper 
et al. 1995), because the latter observations averaged over 
the whole region. 



4.2. Kinematic structures of the molecular gas in 
NGC6334 I 

The striking edge-on orientation between the UCHn re- 
gion and the molecular peaks in NGC6334 1 allows a closer 
analysis of possible kinematic signatures the UCHn re- 
gion might impose on the molecular cores (in addition to 
the heating effect we have described in §4.1.1). Figure 7 
shows two position-velocity (pv) cuts performed for the 
NH 3 (2,2): (a) in east-west direction from the UCHn re- 
gion through the NH 3 (2,2) hyperfine emission hot spot, 
and (b) from the UCHn region in north-western direction 
through the main mm and NH 3 peak position. While there 
is considerable structure in the east- west pv-diagram (a), 
it however does not reveal any clear signature that the 
UCHn could have imposed on the molecular gas. The 
other pv-diagram (b) shows no strong velocity gradients 
from the UCHn region in the direction of the NH3 peak 
position as well. 

Figure 8 presents the moment 2 maps (velocity disper- 
sion maps) of the two NH 3 (1,1) & (2,2) outermost satel- 
lite hyperfine lines and one of the CH3OH lines. While 
the velocity dispersion of the (1,1) hf-line peaks toward 
the north-eastern NH 3 peak, the (2,2) velocity disper- 
sion peaks further south toward the UCHn region and the 
NH3(2,2) hyperfine hot spot. In contrast, the moment 2 
map of CH3OH shows two peaks both associated with the 
two molecular peak positions. 

At estimated temperatures of ~ 100 K the critical densi- 
ties of the CH3OH lines are of the order 10 4 cm -3 , whereas 
the critical densities of the NH3 lines are lower (of order 
10 3 cm~ 3 ). The velocity dispersion of the higher density 
tracer CH 3 OH is dominated by internal motion due to the 
embedded protostars. To some degree this is also the case 
for NH3, and especially the (1,1) moment 2 maps peaks 
toward the strongest mm continuum peak. However, there 
is no obvious (1,1) velocity dispersion peak toward the sec- 
ondary mm peak position. The NH 3 (2,2) velocity disper- 
sion is clearly dominated by the NH 3 (2,2) hyperfine hot 
spot close to the UCHn region. Since NH3 has lower crit- 
ical densities than CH 3 OH, its line emission shows com- 
parably less contribution from the densest cores, and it 
is thus more likely to exhibit temperature and kinematic 
signatures imposed by the adjacent UCHn region. 

4.3. Molecular absorption toward the UCHu region 

In addition to the molecular line emission, there are 
prominent CH 3 OH line absorption features observed in all 
three lines toward the UCHn region in NGC6334I (Fig. 
9). Examining the NH3 data, we find no absorption in 
the satellite hyperfine lines but weak NH3 absorption fea- 
tures in the two main lines. Position- velocity cuts of one 
of the CH3OH and the two NH 3 (1,1) and (2,2) main lines 
through the UCHn region at ADec.=0 reveal that the ab- 
sorption features are blue-shifted with regard to the (weak) 
emission features toward the UCHn region. This signa- 
ture is more pronounced in the CH3OH and NH3(1,1) line 
whereas it is not obviously discernible in the NH3(2,2) line. 
Nevertheless, these data indicate expanding molecular gas 
around the UCHn region. 

5. DISCUSSION 
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The differences between the two massive star-forming 
regions NGC6334I & I(N) in their molecular NH 3 and 
CH3OH emission are rather astonishing: compact thermal 
emission is observed in both species toward the UCH11 re- 
gion/hot core NGC6334I whereas extended NH3 and mas- 
ing CH3OH emission is found toward NGC6334I(N). It is 
not entirely obvious what is causing these differences be- 
cause both sources harbor compact mm continuum sources 
as well as jet-like SiO outflows (Hunter et al. in prep.). 
Since the southern source NGC6334I has a UCH11 region, 
it is likely that evolutionary processes are at least partly 
responsible for the different molecular signatures. How- 
ever, with the data obtained during the present investiga- 
tion we still cannot accurately determine temperature and 
density in both regions, and one expects that the molecu- 
lar emission is sensitive to these parameters. Therefore, it 
is an important next step to get thorough estimates of the 
characteristic temperatures and densities of both regions 
to get a deeper understanding of the physical processes 
causing the observational differences. 



is not the case for NGC6334I(N). There, the gas and 
dust temperatures appear to have decoupled. Addition- 
ally, Sandell (2000) estimated the bolometric luminosities 
of both sources to differ by about two orders of magni- 
tude with 2.6 x 10 5 L Q and 1.9 x 10 3 L Q for NGC6334I 
& I(N), respectively. This difference could either imply 
that NGC6334I(N) may not form a cluster as massive 
as NGC6334I, or NGC6334I(N) may be less evolved and 
younger than NGC6334I and thus may still build up its 
final larger luminosity. With the observations so far, we 
cannot clearly judge which of these two regions is in an 
earlier evolutionary stage. With the UCHn region and 
the infrared cluster present, the whole region NGC6334I 
is older than the region NGC6334I(N), but it is less clear 
whether the same conclusion holds for the mm/NH3 cores 
within NGC6334I. However, it is worth mentioning that 
the non-detection of compact NH 3 emission toward the 
mm sources in NGC6334I(N) is a rarely observed phe- 
nomenon and may indicate a short-lived period during the 
evolution of massive star-forming regions. 



5.1. Four stages of massive star formation within 1 pc 

Figure 11 presents a large-scale overview of the north- 
eastern end of the NGC6334 region comprising NGC6334 1, 
NGC6334E and NGC6334I(N). It is interesting that 
within 1 pc distance we find four different evolutionary 
stages of massive star formation: the shell-like H11 re- 
gion NGC6334E, the champagne- flow- like UCHn region 
NGC6334I and its adjacent outflow driving mm cores, 
and further north the enigmatic source NGC6334I(N). 
NGC6334E lies just outside each of our primary beams, 
thus we have no new data about this source. 

From an evolutionary point of view, the shell-like Hn 
region NGC6334E appears to be the most evolved and 
thus oldest source within the region (e.g., Rodriguez et al. 
1982; Tapia et al. 1996; Carral et al. 2002). The UCHn 
region NGC6334I is significantly smaller than NGC6334E 
and thus apparently in an earlier evolutionary stage (e.g., 
Persi & Ferrari-Toniolo 1982; de Pree et al. 1995). Based 
on the cold dust temperatures, the detection of only a few 
very red infrared sources and weak cm emission (Gezari 
1982; Tapia et al. 1996; Sandell 2000; Carral et al. 2002), 
NGC6334I(N) was always considered to be the youngest 
region within the large-scale NGC6334 complex. However, 
comparing the mm continuum/NH3 sources in NGC6334I 
with those in NGC6334I(N), it is not entirely obvious 
which of these two regions is evolutionary younger. Both 
complexes have compact mm dust continuum sources and 
drive collimated outflows. In contrast to these similari- 
ties, the differences in the compact versus extended NH3 
emission and thermal versus maser CH3OH emission at 
24.9 GHz are striking. Furthermore, the differences in 
the estimated dust temperatures are large with 100 K 
for NGC6334I and only 30 K for NGC6334I(N) (Sandell 
2000). Contrary to this, even constraining only lower lim- 
its for the gas temperatures based on our NH 3 observa- 
tions, we derive kinetic gas temperatures > 100 K for both 
regions. Even considering the different beam sizes of the 
SCUBA observations of the dust emission (between 8.5" 
and 18.5" depending on the frequency) and our NH 3 gas 
observations, the dust and gas temperatures in NGC6334I 
appear to be consistent with each other, whereas this 



5.2. Turbulence versus organized motion 

Figure 12 shows a spectrum of the outermost satel- 
lite hyperfinc component of the NH3(1,1) line in the uv- 
domain on the shortest baseline of 31 m. Several velocity 
components are present: one centered at ~ —9 kms -1 , 
one at ~ — 6kms~ 1 , one at ~ — 3kms -1 , and one at 
~ —1.5 kms . Analyzing the spatial distribution of the 
different velocity components (Figure 13) we find that only 
the — 9 kms -1 component is clearly associated with the 
south-western NH 3 emission peak, whereas all other fea- 
tures are associated with the north-eastern peak (which is 
also the peak of stronger mm continuum emission). 

Contrary to the complex profile in the uv-domain (Fig- 
ure 12), the hypcrfine components of the spectra taken 
individually toward each peak position resemble single- 
peaked Gaussians (Figures 4 & 6). We conducted Gaus- 
sian line-fits to the outermost satellite hyperfine compo- 
nents of the NH 3 (1,1) line toward the two NH 3 peak 
position in NGC6334I and the mm peak position in 
NGC6334I(N), and the derived line- widths Av vary be- 
tween 4.1 and 5.1 kms -1 (see Table 3). 

What is causing these broad line-widths, turbulent mo- 
tion due to, e.g., molecular outflows, or organized mo- 
tion due to, e.g., infall, envelope rotation, disks, or mul- 
tiple components? Sridharan et al. (2002) have shown 
that the NH3 linewidth observed with single-dish tele- 
scopes toward High-Mass Protostellar Objects (HMPOs) 
is significantly lower than toward typical UCHn regions 
and hot cores (Awhmpo ~ 2.1 kms -1 versus Afhotcoro ~ 
3.1 kms -1 ). Since we find gas temperatures of the or- 
der 100 K, NGC6334I & I(N) are considered to be part 
of the hot core group. Therefore, most likely the large 
line-widths are an evolutionary signature of the hot core 
stage and no remnant from the original molecular core. 
The thermal linewidth at 100K is only 0.31 kms -1 and 
thus non-thermal motion has to account for nearly all the 
observed linewidth. To estimate the corresponding en- 
ergy within the core we calculated E = Av 2 t with m 
the core mass derived from the mm con tinuum emission 
(Hunter ct al. in prep.) and Av nt = y/ Av 2 - Aw t 2 hcrmal 
the non-thermal linewidth component. The derived en- 
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ergy values corresponding to internal motion are of the 
order 10 46 erg (see Table 3). 

Interestingly, the typical energy supply of massive 
molecular outflows is of the same order (e.g., Beuther 
et al. 2002). The spatial scales for the molecular out- 
flows arc usually of the order 1 parsec whereas the NH 3 
line- widths observed toward NGC6334I are confined to 
a few thousand AU, thus the estimated energies are not 
directly comparable. However, the uncertainties in deriv- 
ing outflow energies are usually high, up to an order of 
magnitude (e.g., Cabrit & Bertout 1990), and it might 
be possible that comparable amounts of energy could be 
supplied from the outflow on spatial scales as small as 
the NH3 cores. Therefore, turbulence caused by molecu- 
lar outflows may contribute to the observed linewidth in 
NGC6334I. Contrary to this, Zhang et al. (1998, 2002) 
find broad NH 3 line- widths of the order 5kms _1 toward 
two young massive star-forming regions, whereas the NH3 
emission in the surrounding cores has a narrower linewidth 
of the order 2kms _1 . These observations are interpreted 
as indirect evidence for rotating disks within these mas- 
sive star- forming regions. Since we are unable to derive 
any systematic velocity gradient within the separate cores 
of NGC6334I, we cannot attribute our observations to a 
disk scenario. However, it is likely that the observed line- 
widths are caused at least to some degree by organized 
motion within the inner few thousand AU caused by either 
infall, envelope rotation, disks or multiple components, or 
an interplay of various of these components. Only higher 
angular resolution observations can disentangle all these 
contributions. 

5.3. Shocks in NGC6334 1(N) 

In contrast to NGC6334I, the situation appears to be 
different in NGC6334I(N) because the NH3 emission is 
not confined to the mm peak position, and we observe 
broad line-widths over a large spatial area. Since there 
are at least two molecular outflows in the region, the broad 
linewidth could well be caused by shock interactions of the 
outflow with the surrounding molecular core. 

Additional evidence for outflow-shock interactions is the 
detection and the spatial spread of the Class I 25 and 
44 GHz CH3OH maser emission which is believed to be 
caused by the interactions of outflows from massive star- 
forming regions with the surrounding cloud cores (e.g., 
Menten 1991; Sobolev 1993). Class I CH 3 OH maser ob- 
servations at 95 GHz in DR21 revealed that they are lo- 
cated at the interface of the molecular outflow and the sur- 
rounding ambient gas (Plambeck & Menten 1990). Simi- 
lar results were found for Class I CH3OH maser emission 
near 25 GHz in the Orion-KL region, where Johnston et al. 
(1992) suggest that the Class I maser emission is formed 
behind shock fronts between the outflows and the dense 
cloud gas. 

In the case of NGC6334I(N), the north-eastern 25 GHz 
maser feature (the — 3kms _1 component) is obviously spa- 
tially associated with the jet-like SiO(5-4) outflow (Hunter 
et al. in prep., e.g., Fig. 2), and thus fits well into the 
shock-excitement picture. The southern 25 GHz maser 
(the — 5kms _1 component) is less clearly associated with 
one of the two outflow axes. However, the north-west 
south-eastern outflow of Megeath & Ticftrunk (1999) is 



significantly less collimated than the newly detected SiO 
jet of Hunter et al. (in prep.), and therefore it is plau- 
sible that the Class I maser features near the southern 
— 5kms _1 maser peak (i.e., 25 and 44 GHz masers) are 
associated with the interface of the broader outflow cone 
and the ambient cloud. Comparing the overall spatial dis- 
tribution of the 25 and 44 GHz Class I CH 3 OH maser emis- 
sion with the present outflow morphologies (see Fig. 2) all 
maser features appear to be associated with the outflows 
and are likely caused by the shock-interactions of the out- 
flows and the ambient gas. 

The Class I CH 3 OH masers in NGC6334I(N) thus show 
the same behavior as masers of this kind in many other re- 
gions: First, individual narrow-line maser components arc 
spread over areas of several tenths of a parsec, offset from 
any possible exciting source. Additional evidence places 
them in the interaction zones of outflows and dense am- 
bient gas. Second, apart from perhaps Orion-KL, masers 
in the 44 GHz transition are always stronger than those in 
any of the 25 GHz lines, and frequently the strongest of 
all Class I CH3OH masers. This can be understood from 
the (collisional) pumping mechanism of Class I methanol 
masers, which produces efficient masing in the 25 GHz lines 
over just the range of physical parameters one may ex- 
pect in the scenario described above (Leurini, Menten, & 
Schilke 2005, in prep.). This means that possible 25 GHz 
emission from many of the spots from where 44 GHz emis- 
sion arises may be too weak to be detectable. 

6. SUMMARY 

The NH 3 (1,1) & (2,2) and CH 3 OH observations of the 
two massive twin cores NGC6334I & I(N) reveal signifi- 
cantly different emission signatures between both regions. 
While NGC6334I shows strong, compact thermal emission 
in the observed NH 3 and CH 3 OH lines adjacent to the 
UCH11 region and coincident with mm continuum peaks 
observed by Hunter et al. (in prep.), there is no compact 
thermal line emission toward NGC6334I(N). Here, we find 
that the NH 3 emission is distributed over a broad area of 
size approximately 1'. No thermal CH3OH emission is 
found toward NGC6334I(N) but we do detect two Class 
I CH3OH maser peaks with peak brightness temperatures 
of ~ 7000 K. The two maser peaks may be spatially as- 
sociated with 44 GHz Class I CH 3 OH maser emission and 
are likely produced in the interfaces between the molecular 
outflows and the ambient gas. 

NH3 peak brightness temperatures > 70 K toward the 
mm-peak positions in NGC6334I & I(N) indicate kinetic 
temperatures of the same order. However, the NH 3 (1,1) & 
(2,2) lines are only suitable as a good thermometer below 
50 K and thus we cannot derive more than a lower limit for 
the temperatures in these regions. Directly at the western 
edge of the UCH11 region we find a hot spot where the line 
ratios between the NH 3 (2,2) and (1,1) outermost satellite 
hyperfine components are surprisingly high around 4, in- 
dicating even higher temperatures. Since this position is 
not associated with any mm continuum peak, the heating 
is likely caused by the UCH11 region. NH3 emission is also 
detected toward the outflow in NGC6334I, and we can 
estimate a rotational temperature of approximately 19 K 
there. 

Toward the UCH11 region in NGC6334I we observe 
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CH3OH and NH3 absorption. The absorption features are 
blue-shifted with regard to the emission features close-by. 
This signature is consistent with expanding molecular gas 
around the UCH11 region. We do not find additional strong 
kinematic effects imposed from the UCH11 region on the 
adjacent molecular cores. 

We like to thank Tyler Bourke and Paul Ho for help- 
ful ATCA and NH 3 discussion during the analysis and 
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edges financial support by the Emmy-Noether-Program 
of the Deutsche Forschungsgcmeinschaft (DFG, grant 
BE2578/1), and S.T. is grateful to the Alexander von 
Humboldt foundation for a Feodor Lynen research fellow- 
ship. 
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Fig. 1. — Centimeter continuum (grey-scale and contours in the top left image, grey-scale in the other panels), integrated main hyperfine 
component NH3 (1,1) (top right), (2,2) (bottom left) and CH30H(32,i — 3i,2)E integrated line emission toward NGC6334I. The arrows in the 
top-right panel outline the direction of the red and blue outflow lobes (which arc apparently a little bit bent, Hunter et al. in prep.). The contour 
levels are from 10 to 90% (step 10%) from the peak emission of 1197 mjy beam -1 for the cm continuum image (3<r ~ 42 mjy beam -1 ), and in 
3<r steps for the line emission. The approximate 3<r values for the NH3(1,1), (2,2) and CH3OH lines are 21mjybeam~ 1 , 12mjybeam _1 and 
lOmJybeam -1 , respectively. The symbols and spectral lines are labeled in each panel. The synthesized beams are shown at the bottom-right 
of each panel as well. 
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Fig. 2. — Integrated main hyperfine component NH3 (1,1) and (2,2) line emission toward NGC6334I(N) (left and right, respectively). 
The two lines outline the directions of the two observed outflows in the region (north-east south west Hunter et al. in prep., south-east 
north west Megeath & Tieftrunk 1999). The contour levels are 58(58)348 mjy beam -1 (3a ~ 9mjybeam _1 ) and 33(33)198 mjybcam -1 
(3cr ~ 6mjybeam _1 ) for the NH3(1,1) and (2,2), respectively. The symbols are marked at the bottom-right of each panel. The synthesized 
beams are shown at the bottom-left of each panel as well. 
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Fig. 3. — CH3OH maser spectra taken toward the NGC6334I(N). The top-left spectrum shows the uv-spectrum on the shortest baseline of 
31 m presenting all detected components. The top-right spectrum is taken after imaging toward the north-eastern CH3OH maser position (the 
conversion factor for K/(Jy beam -1 ) of the right spectrum is ~ 385). The two bottom spectra are the same uv-spectra as the top-left one but 
this time presented on the velocity-scale and differentiated between the CH30H(32,i — 3i ; 2)E and 0113011(42,2 — 4i 3)E lines (bottom-left 
and bottom-right), respectively. 
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Fig. 4. — NH3 spectra toward the two peak positions in NGC6334I (the offsets of the phase reference centers listed in Table 1 are marked 
in each panel). The rather poor fits are shown in dotted lines. 
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Fig. 5. — Integrated images of the outer satellite hypcrfinc components at the most negative velocities toward NGC6334I. The left and 
middle panels show satellite hyperfine components of Nfi3(l,l) (the blend of the NH3(1,1).F, F\ = 0.5, — 0.5, 1 & F, Fi = 0.5, — 1.5, 1 lines 
at ~ 23.696 GHz) and (2,2) (the NH 3 (2,2)F = 2-1 line at ~ 23.725 GHz)), respectively. The contour levels for the NH 3 (1,1) & (2,2) satellite 
components are in 3<x steps with 3<r(l, 1) ~ lOmJybcam -1 and 3<r(2, 2) ~ 9mjybeam~ 1 . The right panel presents a ratio map of the two 
integrated satellite hyperfine images (2,2)/(l,l). The cm continuum emission is presented in grey-scale in the left and middle panels and in 
contours (10 to 90% from the peak emission of 1197 mjy beam - x ) in the right panel. The symbols are the same as in Figure 1. 
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Fig. 6. — NH3(1,1) and (2,2) spectra taken toward the nominal position of the mm continuum peak discovered by Hunter et al. (in prep.) 
toward NGC6334I(N). The offsets of the phase reference centers listed in Table 1 are marked in each panel. The fits arc shown in dotted 
lines. 



14 



Bcuther et al. 



-12-000 - -5.500 




3 


c 












velocity (km/s) 



Relative R.A. foresee) 

Fig. 7. — Position- velocity cuts through the NH3(2,2) outer satellite hyperfine line data cube (NH3(2,2)_F = 2 — 1, integrated from —12 to 
— 5kms -1 ) of NGC6334I along the axes outlined in the left panel. The middle panel corresponds to the cut in east-west direction, the right 
panel shows the cut north-west south-east direction. The contour levels are in 11.5, 5.0 and 16.5 mjy beam -1 steps for the left, middle and 
right panel, respectively. 
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Fig. 8. — The grey-scale shows moment 2 maps of the NH3(1,1) (left), (2,2) (middle) outermost satellite hypcrfine lines and the CH30H(32,i — 
3i.2)E line (right panel) in NGC6334I (the wedge is in units of kms -1 ). The contours in all three panels show the cm continuum data from 
10 to 90% of the peak emission of 1 197 mjy beam -1 . The symbols are the same as shown in Fig. 1. 
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Fig. 9. — CH3OH spectra taken toward the two emission peaks and the UCH11 region (right panel) in NGC6334I. The offsets are given 
with respect to the phase center listed in Tabic 1. 
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Fig. 10. — Position-velocity cuts in the R.A. direction at ADec.= 0" for the CH 3 OH(32,i - 3i,2)E line (left) and the two NH 3 (1,1) and 
(2,2) main components (middle and right) in NGC6334I. Positive emission is in full contours and negative emission in dashed contours. The 
contour levels are in ±10, ±10 and ±9.5 mjy beam -1 steps for the left, middle and right panel, respectively. 
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Fig. 11. — Large-scale overview of the northern end of the NGC6334 complex. The grey-scale shows the cm continuum emission of 
NGC6334I and the contours present the integrated NH3(1,1) emission of the main hyperfine component. The contour levels are from 10 to 
90% of the peak values of each image with peak values: cm emission 1197 mjy beam" 1 , NH 3 (1,1) in NGC6334I 175 mjy beam" 1 , NH 3 (1,1) 
in NGC6334I(N) 386 mjy beam -1 . The grey circle and star present the approximate structure of NGC6334E (e.g., Carral et al. 2002). 
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Table 1 
Phase reference centers 



Source 


R.A. 


Dec. 




[J2000] 


[J2000] 



NGC6334I 17:20:53.44 -35:47:02.2 
NGC6334I(N) 17:20:54.63 -35:45:08.9 



Table 2 

CH3OH Class I maser positions 



Velocity 


R.A. 


Dec. 


[kms-i] 


[J2000] 


[J2000] 


-3 


17:20:55.68 


-35:45:00.2 


-5 


17:20:55.03 


-35:45:15.0 



Note. — The estimated absolute po- 
sitional accuracy is <~ 0.2". 



Table 3 

NH3(1,1) satellite hyperfine components 



Source 


Offset 


Av 


E 




H 


[kins- 1 ] 


[erg] 



NGC6334I -1.0,4.5 5.1 1.2 x 10 

NGC6334I -3.0,3.5 4.1 4.5 x 10 

NGC6334I(N) 6.5,5.0 4.5 4.8 x 10 



Note. — The line-widths are derived from Gaus- 
sian fits to the outermost NH 3 (1,1) hyperfine com- 
ponent at most negative velocities (see Figures 4 & 
6). 



